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We present polaron effective masses and selected polaron band structures of the Holstein molec- 
ular crystal model in 1-D as computed by the Global-Local variational method over a wide range 
of parameters. These results are augmented and supported by leading orders of both weak- and 
strong-coupling perturbation theory. The description of the polaron effective mass and polaron band 
distortion that emerges from this work is comprehensive, spanning weak, intermediate, and strong 
electron-phonon coupling, and non-adiabatic, weakly adiabatic, and strongly adiabatic regimes. Us- 
D ' ing the effective mass as the primary criterion, the self-trapping transition is precisely defined and 

5^ , located. Using related band-shape criteria at the Brillouin zone edge, the onset of band narrowing 

is also precisely defined and located. These two lines divide the polaron parameter space into three 
regimes of distinct polaron structure, essentially constituting a polaron phase diagram. Though the 
self-trapping transition is thusly shown to be a broad and smooth phenomenon at finite parameter 
values, consistency with notion of self-trapping as a critical phenomenon in the adiabatic limit is 
demonstrated. Generalizations to higher dimensions are considered, and resolutions of apparent 
'""O ■ conflicts with well-known expectations of adiabatic theory are suggested. 

O ' PACS numbers: 71.38.+i, 71.15.-m, 71.35.Aa, 72.90.+y 
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I. INTRODUCTION 

m _ 

The Holstein model is experiencing a resurgence of in- 
terest due in part to progress in the modelling of high- 
T c oxides and to advances in techniques of theoreti- 
cal analysis, including variational methods llHfj], density 
Q\ , matrix renormalization group (DMRG) , quantum 
^ ■ Monte Carlo (QMC) |P-[l5[|, cluster theor y |16[-p4[| , weak- 
coupling perturbation theory (WCPT ) |25H30[| , strong- 
coupling perturbation theory (SCPT) |ggg|l|-|||. 

Our principal quantitative tool is the Global-Local 
variational method, supported at weak and strong cou- 
pling by the leading orders of perturbation theory. Else- 
where we have made a number of specific quanti- 
tative comparisons with other high-quality methods at 
general points in the polaron parameter space, demon- 
strating the quantitative consistency of our results with 
the best available by any method, and with all known 
limiting behaviors, whether at strong or weak coupling, 
or in the adiabatic or non-adiabatic limits. 

There is an important exception to this broad consis- 
tency among independent methods. A number of well- 
known results, important in both their quantitative and 
qualitative aspects, are associated with the adiabatic ap- 
proximation J34]-[49||. In particular, the concept of the 
self-trapping transition that arises most universally from 
adiabatic theory is that of an abrupt change from infinite- 
radius or "free" states at weak coupling to finite-radius 
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or "self-trapped" states at strong coupling. Within the 
adiabatic perspective it is thus asserted that there is no 
self-trapping transition in 1-D (only in 2-D and 3-D) since 
adiabatic theory generally finds 1-D polaron states to 
be characterized by finite radii at any coupling strength. 
Although these finite-radius states may be described as 
either "large" polarons or "small" polarons depending 
on coupling regime, no essential distinction is recognized 
since polaron structure is found to be essentially the same 
for all coupling strengths. 

This adiabatic perspective on self-trapping is not sup- 
ported by our results. By examinging polaron proper- 
ties with an accuracy and scope previously unavailable, 
we are able to demonstrate that despite a partial consis- 
tency with adiabatic theory, a partial inconsistency with 
adiabatic theory persists up to and including the adi- 
abatic limit. In particular, we find a precisely-definable 
self-trapping transition to exist in 1-D separating polaron 
states of distinct structure at strong and weak coupling, 
and we find that these distinctions persist into the adia- 
batic limit. Our demonstration of these findings here is 
limited to those aspects that are discernible from the po- 
laron effective mass and overall band distortion; however, 
more detailed and fully-corroborating findings based on 
analyses of internal polaron structure are presented else- 
where [HJ3^]. 

Our purpose here is to examine a large volume of 1- 
D variational results as empirical data of high but less 
than ultimate precision, and to subject that data to the 
synthetic exercise of formulating globally-consistent in- 
ferences about the exact underlying physics of the prob- 
lem. Some of these inferences can be tested against other 
approximate results or known limits, but the most in- 
teresting perhaps are those for which there is not nec- 
essarily any formal demonstration available. Near the 
close of our paper, we augment these 1-D results with in- 
dependently ascertainable information regarding dimen- 
sional relationships, lifting certain of these inferences into 
higher dimensions. 

We are particularly interested in quantifying charac- 
teristics of polaron band shape, including, for example, 
the polaron effective mass and polaron band width. Us- 
ing such, we identify objective criteria that permit the 
accurate delineation of phase boundaries that together 
constitute the polaron phase diagram. The essential fea- 
tures of this phase diagram are: 

i) the self-trapping line, gsr, separating the small po- 
laron regime from the intermediate and large polaron 
regimes. 

ii) the line indicating the onset of band narrowing, g^, 
separating the large polaron regime from the intermedi- 
ate and small polaron regimes. 

iii) the termination of the self-trapping line at finite 
coupling in the non-adiabatic limit. 

iv) the termination of the onset line at the non- 
adiabatic/adiabatic crossover in the weak-coupling limit. 

v) the convergence of the self-trapping and onset lines 
at the self-trapping critical point in the adiabatic limit. 
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While the adiabatic critical point has been known for 
some time, none of the other above-noted features of the 
polaron phase diagram have been determined with quan- 
titative precision, and some of these features have not 
existed as theoretical constructs prior to our present anal- 
ysis. 

This paper is organized as follows: In Section II, we 
present the model and states upon which the present 
work is based, and set down notation. In Section III, 
we focus on the global ground state energy, displaying 
specific results according to our own method and com- 
paring our results with those of other authors and certain 
approximate formulas. In Section IV, we turn to the po- 
laron effective mass, computing effective mass curves and 
analyzing these to determine the precise location of the 
self-trapping line. In Section V, we present a number 
of typical polaron energy bands and perform a detailed 
study of the dependence of the polaron bandwidth on 
system parameters, permitting a characteristic line to be 
determined that marks the onset of polaron band nar- 
rowing. In Section VI, we synthesize our results in the 
form of a phase diagram encompassing the entire prob- 
lem, characterizing distinct polaron "phases", the "tran- 
sitions" between them, and extracting the critical behav- 
ior in the adiabatic limit. In Section VII we discuss the 
implications of our results to polaron problems in higher 
real-space dimensions. Conclusions are summarized in 
Section VIII. 



II. MODEL AND STATES 

As our system Hamiltonian, we choose the traditional 
Holstein Hamiltonian f36|j5Ci|| 

H = H ex + H ph + H ex ~ vh , (1) 
H ex = E^2aia n - J^at(a„+i + a n _i) , (2) 

n n 

HP h = fuvJ2blb n , (3) 



H ex ~ ph = -ghoj aW(bl + b n ) , (4) 

n 

in which creates an exciton in the rigid-lattice Wannier 
state at site n, and b^ creates a quantum of vibrational 
energy in the Einstein oscillator at site n. We presume 
periodic boundary conditions on a one-dimensional lat- 
tice of N sites. The exciton transfer integral between 
nearest-neighbor sites is denoted by J, uj is the Einstein 
frequency, and g is the local coupling strength. (Except 
where displayed for clarity, the reference energy E is set 
to zero throughout.) These are dimensioned quantities, 
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from which two dimensionless control parameters can be 
defined in different ways. 

One non-dimensionalizing scheme involves selecting 
the phonon quantum Tim as the unit of energy; in these 
terms, the natural dimensionless parameters are the elec- 
tronic hopping integral (J/Tnu) and the exciton-phonon 
coupling constant (g). This scheme is particularly appro- 
priate when considering dependences on J and/or g at 
fixed uj, such as we shall be concerned with in most of 
this paper. 

Near the adiabatic limit, many common considerations 
involve both J/Hlo and g as being large or diverging in a 
certain asymptotic relationship, for which purposes it is 
sometimes convenient to non-dimensionalize by selecting 
zJ as the unit of energy, where z is the site coordina- 
tion number (z — 2 in 1-D). In these terms, the nat- 
ural dimensionless parameters are the phonon quantum 
(7 = hu)/2J) and the "small polaron binding energy" 
A = e p /2J = g 2 hui/2J. The adiabatic limit is reached 
by allowing 7 to vanish at arbitrarily fixed A. We dis- 
tinguish these two conventions as non-adiabatic (unit = 
hcu) and adiabatic (unit = zJ) scaling, respectively; ex- 
cept where explicitly noted, we conform to non-adiabatic 
scaling. These conventions are not to be confused with 
other, regime-specific terms; for example, throughout 
this paper, we use the term non-adiabatic regime to mean 
J/huj < 1/4, adiabatic regime to mean J/Tiuj > 1/4, and 
adiabatic limit to refer to an extreme limit J/Tllo — ► 00. 

Our central interest in this paper is in the polaron en- 
ergy band, computed as 

E{k) = (9(k)\H\*(k)) , (5) 

wherein H is the total system Hamiltonian, |^(/t)) is 
a normalized trial Bloch state, and n is the total joint 
crystal momentum label of the exciton-phonon system. 
The set of E(k) so produced constitute an estimate (up- 
per bound) for the polaron energy band |5^j5^]. All of 
our calculations are performed using the variational trial 
states of the Global- Local method 0: 

|*(«)) = |«)/<«| K ) 1/2 , (6) 



|«) = ]T e to <-»4 exp{-iV-3 

nn a \_ q 

[(^e-^-^e^)bl~H.c.}}\0) , (7) 

such that polaron structure is represented through three 
classes of variational parameters {a%, P*,^}. These 
states are eigenfunctions of the appropriate total momen- 
tum operator and orthogonal for distinct k with the re- 
sult that H is block-diagonal in k, making variations for 
distinct k independent plj . The self-consistency equa- 
tions that follow, the method of solving those equations, 
and sample results have been detailed in ||. 

In our approach, a "complete" variational solution con- 
sists of a set of N variational energies E(k) and N po- 
laron Bloch states \^(k)), the latter being described by 
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a distinct set of variational parameters for each k. In 
all, our original data in this paper are drawn from nearly 
1200 complete band structure calculations, though for the 
most part only selected energies or derivatives thereof 
arc displayed in large graphical compilations. These 
data roughly span a rectangle of the (J/Ttui, g) parameter 
space from the origin out to the extreme point (9,4.5). 
Though both small- J /%uj and small-g calculations were 
performed, our main interest has been in J/Hlo > 1/4, 
and g > 1/2. In terms of the adiabatic parameters 7 and 
A, the figure covered by our data appears to be substan- 
tially larger; the parameter 7 quantifying adiabaticity 
ranges over two orders of magnitude (0.056, 5), and the 
parameter A quantifying coupling strength ranges over 
nearly five orders of magnitude (0.00056, 40.5). The data 
presented herein thus include essentially every regime 
subject to analytical and numerical analysis, with the 
technical exception of certain limits beyond the reach 
of finite-parameter sampling, implying a scope that is 
greater than has been assessed heretofore by any tech- 
nique of comparable reliability. 

III. GROUND STATE ENERGY 

Our approach divides the total Hilbert space of the 
problem into subspaces of the total crystal momentum 
within which the ground state of each k sector is deter- 
mined independently. In this sense, the global ground 
state energy, E(0), is only one of N independent ground 
state energies computed on an equal footing. There is no 
guarantee that the numerical value of E(0) constitutes 
any more (or less) accurate an estimate of its particular 
target value than any other E(n) so determined. (Anec- 
dotal evidence from convergence characteristics does sug- 
gest that there is some K-dependence to numerical accu- 
racy in our computations; however, it is not true that 
k = is always the best-converged k value.) 

Figure [l] shows the dependence of the global ground 
state energy E(0) on the coupling strength for various 
values of J/hw. The overall behavior of the ground state 
energy is to trend between two asymptotic g 2 depen- 
dences with differing coefficients and offsets. 
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FIG. 1. The global ground state energy E(0)/hui vs. 
g 2 . Solid curves are the Global-Local variation energies for 
J/hcu = 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, and 9.0. The diag- 
onal line is — g 2 , the limiting value of the polaron energy as 
J/Tiu — » for any g and as g — » oo for any J. The straight 
dotted lines asymptoting each Global-Local curve at small g 2 
are given by the WCPT formula ([]) . The arched dashed lines 
asymptoting each Global-Local curve at large g 2 are given 
by the SCPT formula Diamond symbols indicate the 

self-trapping points as determined in Section IV. Square sym- 
bols indicate the onset of band narrowing as determined in 
Section V. 



Using weak-coupling perturbation theory (WCPT), 
one can show that the leading dependence of the ground 
state energy on the coupling constant is given for any 
J /hut by 

S(0)«-2J- , ^ . (8) 

This is, in fact, what we find to within numerical preci- 
sion within the weak coupling regime of the Global-Local 
method. 

On the other hand, the strong-coupling perturbation 
result (in one dimension) 

E(0)*-2Ja(i + ±)=- 9 *Hu J -^- (9) 

that provides the correct asymptotic behavior at large 
A does not offer any improvement over the WCPT esti- 
mate until g gsT- This result breaks down definitively 
for A < 1/2, where its dependence on the coupling con- 
stant fails to capture even the correct qualitative trend. 
(Except for very small J/hu), this characteristic is not 
significantly changed by using the complete second-order 
expression shown in ( |12|) below.) 

The "knee" between weak- and strong-coupling trends 
in E(0) bears the hallmarks of the self-trapping tran- 
sition and could be used (as we have done in Ref. ||) 
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as a locating criterion. Here, however, we focus more 
broadly on polaron band shape and how the character- 
istic fc-dependent distortions of the band are organized 
relative to the fundamental self-trapping event. On fig- 
ures throughout this paper, as in Figure 0, we have indi- 
cated two significant sets of self-trapping-related values 
as determined by band-shape criteria in the sections that 
follow. The diamond symbols indicate our location of the 
self-trapping line gsr using the polaron effective mass 
(Section IV), and the square symbols indicate our loca- 
tion of the onset of band narrowing g^ using the band 
edge curvature (Section V). It is clear that the gsr points 
coincide well with the locations of the "knee" in the 
ground state energy, confirming that the effective mass 
and the ground state energy are mutually-consistent lo- 
cators of the self-trapping transition. This is no surprise, 
of course, in view of the intimate relationship between 
the ground state energy and the effective mass. Though 
no dramatic feature exists in the ground state energy at 
the onset points g??, they do appear to correlate with the 
beginning of a turnover from the limiting weak-coupling 
trend into a more intermediate behavior; this proves to 
be generally true, as will be developed in the following. 

IV. EFFECTIVE MASS AND SELF-TRAPPING 

The polaron effective mass is the traditional indicator 
of the self-trapping transition, with the transition most 
commonly being taken to be evidenced by a jump in the 
effective mass upon self-trapping. Although this charac- 
terization is convenient and widespread, it is generally 
expected that this rise in the effective mass should be 
rapid but continuous at finite J/huj and g |53|-|56||. 

The location of the transition has not been well char- 
acterized, in part because this blend of smoothness and 
abruptness tends not to be well captured by most of the 
approximate methods that for many years have been ap- 
plied to the problem. It is common for some approxi- 
mate methods to miss the transition completely, passing 
smoothly into invalidity as the transition region is ap- 
proached and crossed from one or the other side; other 
methods may exhibit signs of breakdown near the tran- 
sition. It is thus that discontinuities, crossings of solu- 
tion branches, intersections of weak and strong coupling 
asymptotics, and other symptoms ultimately but impre- 
cisely related to self-trapping effects have been used as 
rough locators of the fundamental self-trapping event. 
While each of these estimations has its own merit, all 
suffer from the fact that they indicate more directly an 
insufficiency of knowledge than they do the positive an- 
swer actually being sought. Each approach further suf- 
fers from quantitative ambiguity, as is evidenced by the 
wide variety of order-of-magnitude characterizations to 
be found in the literature. 

Since we have achieved substantial improvement in 
the quality of variational results over a wide range of 
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parameters (weak to strong coupling, non-adiabatic to 
adiabatic), we are in position to attempt to locate the 
physically-meaningful self-trapping transition by posi- 
tive, objective criteria applied to quantitative trends of 
measurable physical properties. 

In this regard, we are particularly interested in using 
our methods to determine as precisely as possible the 
behavior of the polaron effective mass. We compute the 
effective mass using the definition 

to* _ 2J 

U Ok 1 l"=° 

using a discrete representation of the k derivative at the 
Brillouin zone center. (Note: the bare effective mass too 
was also computed as a discrete derivative on a finite 
lattice rather than in the N — > oo limit; this is formally 
necessary so that the ratio to*/too on a finite lattice tends 
exactly toward unity in the weak-coupling limit.) 

Looking at the dependence the effective mass on J/Hlu 
and g as computed under the Global-Local method (ac- 
tually, log(m*/mo), see Figure ||), it is clear that for most 
J/hu there are two distinctly different trends in the de- 
pendence of the effective mass on the exciton-phonon 
coupling constant, one at weak coupling the other at 
strong coupling. 
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FIG. 2. The ratio of the polaron effective mass to the bare 
effective mass (m* /mo) is plotted against exciton-phonon 
coupling in a logarithmic scale. Solid curves, left to right: 
J/huj = 0.25,0.5,1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, and 9.0. 
Dashed curve: The exact J/hu = + limit. Open diamonds: 
Self-trapping points as determined by inflection criteria as in 
Fig. b, c, and d. Shaded region: Area lying between up- 
per and lower bounds for self-trapping points as determined 
from distortion characteristics as in Fig. ^ and f. Solid dia- 
monds: Self-trapping points as determined by kinetic energy 
analysis as in Ref. M. Squares: Onset of narrowing points as 
determined by band edge criteria as in Section V. 

Further, these asymptotic dependences are separated 
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by a zone through which a relatively rapid transition oc- 
curs. It is this relatively rapid transition between asymp- 
totic dependences that we identify as the physically- 
meaningful self-trapping transition. We identify the self- 
trapping transition with the point of most rapid increase 
in the effective mass in this transition region, uniquely 
identified mathematically by an inflection point in the in- 
dependence of log(m*/mo) a t fixed J/Tilj. Thus, where 
the numerical data was sufficiently smooth, we applied 
the criterion 



9 st 3 -7-0 l°g I 

dg z \m 



= (11) 

g=gsT 



as indicated by the zero-crossings of the curves shown in 
Figure I for J/Tioj = 2.0, 3.0, 4.0, and 5.0. 
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FIG. 3. Curvature of the logarithm of the effective mass with respect to the exciton-phonon coupling constant in the vicinity 
of the self-trapping transition. J/hw = 2.0, 3.0, 4.0, 5.0, 6.0, and 7.0. 



For J/fnu — 1/4, 1/2, and 1.0 (curvature assays omit- 
ted), the variational band energies were generally well- 
converged and smooth for most purposes; however, no 
inflection feature could be resolved within the limits of 
our computation below J/Huj « 1.5. One can show using 
(HI) ( see Section V) that no inflection feature exists in 
the small J/hu limit; thus, we expect the inflection fea- 
ture to disappear below some finite J/Huj. It is possible 
that this disappearance occurs between J/htv = 1.0 and 
1.5, but we are unable to make that determination on the 
basis of our present data. 

For J/Tiui = 6.0, 7.0, and 9.0, the accurate location of 
the self-trapping transition is hindered by the emergence 
of method-dependent distortions of the transition fea- 
ture. These distortions grow out of the transition point 
itself and develop into full discontinities with increasing 
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adiabaticity as can be seen in the J/Tilo — 7.0 case. The 
J/fruj = 6.0 case falls just below the appearance of the 
first discontinuity. 

That the J/Tuc = 7.0 case exhibits two discontinuities 
is a finite-size effect arising because the band energies at 
K = and k = ±2ir/N (values used to compute the effec- 
tive mass) experience discontinuities at slightly different 
values of g. In an infinite lattice, these values become 
infinitesimally separated such that the two apparent dis- 
continuities combine into one. 

Rather than place inordinate weight on the fine struc- 
ture of artifacts, we interpret self-trapping artifacts as 
rather roughly locating the self-trapping transition and 
use the obvious distortions and discontinuities to define 
generous error bars within which the transition can be 
safely assured to lie. Thus, in Figure we have indi- 
cated effective mass values within these error bars by a 
shaded region. These uncertainties are apparent in the 
effective mass because of the extreme steepness of the ef- 
fective mass at the self-trapping transition; elsewhere in 
this paper, the " error bars" associated with self-trapping 
artifacts are significantly smaller than the plot symbol 
size and so are omitted without further comment. 

The solid-diamond symbols included in the shaded re- 
gion correspond to self-trapping points as determined in 
Ref. (|] using the kinetic energy rather than the effective 
mass as the diagnostic criterion. 

The values of gsr, including estimated errors at higher 
J/hu and values determined from the kinetic energy 
analysis are listed in Table |. 

TABLE I. Values of gsr and qm as determined by the 
Global-Local method. Ranges (•••,•••) indicate estimated 
upper and lower bounds as discussed in the text. * = qst 
values determined by analysis of the kinetic energy in Ref. fe) . 
f = <7jv limiting value inferred from formal considerations (cf. 
(Ell) ff.). £ = gjv value adjusted to compensate for distortion. 
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1.5000 
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2.0000 




1.21562 


2.0 


2.4142 


2.42027 


1.71000 


3.0 


2.7321 


2.71145 


2.0826 


4.0 


3.0000 


2.96890 


2.37364 


5.0 


3.2361 


3.20259 


2.61851 


6.0 


3.4495 


3.422* G (3.416,3.424) 


2.84016 


7.0 


3.6458 


3.6288* G (3.6230,3.6450) 


3.04156 


9.0 


4.0000 


4.012* G (4.010,4.018) 


3.4278} 
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Our effective mass analysis is thus consistent with prior 
estimations of the self-trapping line using the kinetic en- 
ergy, ground state energy, phonon energy, and exciton- 
phonon interaction energy set forth in previous work 

in. 

In specific cases, these values can be compared with 
results of other high-quality methods that have analyzed 
the self-trapping transition as we describe it. In quantum 
Monte Carlo studies, for example, the "break" observed 
in the electron-phonon correlation function for the 1-D, 
J/hu — 1 case coincides well with the value tabulated 
here (see Figure 10 of Ref. ||, with appropriate rescal- 
ing of parameters). Similarly, the break observed in the 
kinetic energy for the same case, though less sharp, is 
also consistent with our placement of the self-trapping 
line (see Figure 4 of Ref. ||). 

V. ENERGY BANDS AND BAND NARROWING 

If there is any one archetypical polaronic effect, it is 
perhaps the strong narrowing of the quasiparticle energy 
bandwidth (J — > J — Je~ 9 ) as reflected in the first order 
of strong-coupling perturbation theory (SCPT); to this 
order and in this narrow view, the polaron bandwidth 
(4Je~ 9 ), effective mass (m*/mo = e 9 ), and ground 
state energy (—g 2 hui — 2Je~ 9 ) are so simply related as 
to appear interchangeable and redundant as characteri- 
zations of the polaron structure. 

This convenient and intuitive picture is quite incom- 
plete, however, in numerous ways. Consider the contin- 
uation of SCPT to second order, which yields ]33| 

E(k) — —g 2 huj 

2 

—2Je 9 cosk 

r [/(2 5 2 ) + /(. 9 2 )cos2 K ] , (12) 

TlbJ 

where f(x) — Ei(a;) — 7 — ha(a^), in which Ei(x) is the 
exponential integral and 7 is the Euler constant. The 
second-order correction makes a strong contribution to 
the ground state energy (second order typically substan- 
tially larger than first order), a weak contribution to the 
effective mass (second order larger or smaller than first 
order), and makes no contribution to the polaron band- 
width (first order uncorrected by second order). The 
widely held, rough picture of polaronic trends and re- 
lationships contained in the first order of SCPT is thus 
quite limited, though accurate at very small J/Jiui and 
useful in a heuristic sense. 

The second order of SCPT as displayed above is also 
quite limited at finite n. Just as the typically-important 
correction to the ground state energy comes not from first 
but second order, one can plausibly infer from the struc- 
ture of SCPT that important corrections to the effective 
mass, bandwidth, and other finite-K properties remain to 
be extracted from third and higher orders of SCPT. This 
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is suggested as well by quantitative comparisons of the 
second-order SCPT effective masses with Global-Local 
data and density matrix renormalization group data as 
presented in Ref. [Q. Correct, however, is the qualitative 
characteristic that the second-order correction modifies 
the polaron band shape, flattening the band somewhat 
in the outer Brillouin zone. 

More generally, away from the strong coupling limit, 
polaron energy bands are more strongly distorted figures 
whose non-sinusoidal dependence on k is a crucial reflec- 
tion of polaron structure. For such bands, the polaron 
binding energy, effective mass, and polaron bandwidth 
no longer stand in any simple relationship to each other. 
In particular, the general independence of the polaron ef- 
fective mass and the polaron bandwidth away from the 
strong coupling limit will be exploited to significant ad- 
vantage in the following. 

In the limit g — > + , one finds the polaron energy band 
assuming a clipped form, 

E{k) = -2Jcos(k) \k\ < k c , (13) 

= -2J + huj \k\ > k c , (14) 

reflecting the difference in the character of polaron states 
above and below the wave vector k c (given by the con- 
dition 2J[1 — cos(kc)] = Tiuo) at which the bare exci- 
ton energy band crosses into the one-phonon continuum; 
when J/Ticu < 1/4, n c does not exist. Although this 
crisply clipped band form is strictly valid only in the limit 
g — > + , its essential characteristics persist to non-trivial 
values of the exciton-phonon coupling strength. Using 
WCPT, for example, one can show that 



E(k) = -2Jcos(k) 

g 2 huj 



y/[hu + 2Jcos(K)} 2 -4J 2 



\n\ < k c . (15) 



This correction diverges for any fixed g as k — > k c , sig- 
nalling the breakdown of this perturbation theory in the 
outer Brillouin zone where a qualitatively distinct behav- 
ior is found. 

Our variational approach does not suffer this problem, 
allowing us to compute polaron energy bands as illus- 
trated in Figure |^ for J/hu> ~ 4 for assorted coupling 
strengths, beginning at the weak coupling limit g = + 
and proceeding through the self-trapping transition into 
the strong coupling regime. 
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FIG. 4. Shifted energy bands [E(k) - E(0)]/hu) for 
J/huj = 4.0, g = 0.4, 0.8, 1.2, 1.6, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 
3.2 as computed by the Global-Local method. k c ~ 0.16097T. 
(Note: This illustration approximates N = 64 by quadrati- 
cally interpolating between N — 32 data.) 

As the exciton-phonon coupling is first turned on, its 
first effect on the energy band is a rounding of the an- 
gular feature at re = ±re c . This rounding can be un- 
derstood in terms of quantum-mechanical level repulsion 
in the vicinity of the crossing of the bare exciton band 
and the edge of the one-phonon continuum. For a sig- 
nificant interval of coupling strength, the primary effect 
of increasing coupling is a progressive rounding of this 
feature and a gradual smoothing of the overall energy 
band that spreads outward from re c to higher and lower re. 
Throughout this weak-coupling phase of energy-band dis- 
tortion, the shape of the polaron energy band is decidedly 
non-cosine, and the polaron bandwidth bears little rela- 
tion to the polaron effective mass. Only after reaching a 
significantly finite value of the exciton-phonon coupling 
strength (g w g^ ~ 2.4) does this smoothing of the en- 
ergy band progress to the point that the overall polaron 
energy band begins to narrow significantly. Once the 
self-trapping transition is encountered (g « gsr ~ 3.0), 
the band is nearly cosine-shaped, approaching perfect co- 
sine form asymptotically at large coupling; in this self- 
trapped regime, the familiar direct relationship between 
the strong-coupling polaron bandwidth and the polaron 
effective mass is recovered. 

A curvature analysis (not shown) identifying the point 
of maximum negative curvature of E(k) with respect to 
re shows that throughout this process, up to the self- 
trapping transition, the "knee" in the energy band re- 
mains located at re » re c . These features indicate rather 
strongly that the notion of distinct outer band (|re| > re c ) 
and inner band (|re| < re c ) structure persists from g = + 
up to the vicinity of the self-trapping transition. 
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A. Onset of Narrowing 



Much as the polaron effective mass is a pure zone- 
center quantity with which to characterize the self- 
trapping transition, we now consider a pure zone-edge 
quantity with which to characterize the onset of band 
narrowing. 

When exciton-phonon coupling is small, and in par- 
ticular in the limit g — > + , the polaron energy band is 
flat at k — ±7r. When exciton-phonon coupling is large, 
and in particular in the limit g — > oo, the polaron en- 
ergy band is again flat at n = ±tt. At intermediate cou- 
pling strengths, the energy band has a nontrivial negative 
curvature at the Brillouin zone edge, implying that this 
curvature passes through an extremum at finite g. In- 
spection of Figure ^ renders plausible the inference that 
such dips in the band edge curvature may be meaning- 
fully associated with the transition from the smoothing 
to the narrowing regime. We have analyzed the curva- 
ture of the polaron energy band at the Brillouin zone 
edge, identifying the point of extreme curvature by the 
relation 

d_ d?_ 

N dg dn 2 

at fixed J/hio. The band-edge curvatures of a number 
of polaron energy bands spanning a large range of J/hco 
and g are displayed in Figure ^. 



(16) 




FIG. 5. Curvature of the polaron energy band with re- 
spect to the exciton-phonon coupling constant evaluated at 
the Brillouin zone boundary. The dashed line at lower left 
is neither formally derived nor computed, but represents an 
extrapolation of our J/Tioj — 1/4 data to the weak-coupling 
limit. The unit of energy is tuo. 

For most values of J/huj, the curvature reversal that 
marks the onset of narrowing is broad in g and easily 
located; the numerical value of g^ can be determined 
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rather precisely by fitting a parabolic curve to the few 
points near the minimum. The values so determined have 
been tabulated in Table |, and are indicated in figures 
throughout this paper by square symbols. This exercise 
demonstrates that the extrema in the band edge cur- 
vature are very effective locators for the onset of band 
narrowing, just as the self-trapping points are very ef- 
fective locators of the culmination of band narrowing. 
There is no outstanding feature of zone-center quantities 
such as the ground state energy or effective mass at the 
onset points, as might be expected since the onset phe- 
nomenon is mostly a band edge effect, although the crite- 
rion g ~ gN does appear at least roughly to characterize 
the coupling scale at which both the ground state energy 
and the effective mass begin to deviate significantly from 
WCPT. 

At J/hu = 9, the curvature reversal data appear to be 
distorted relative to the nearly-symmetric figures found 
at lower J/Tkv. We interpret this as a method-dependent 
artifact, not unlike nor unrelated to the jump discontinu- 
ity in the effective mass also found at this J/fiw. For this 
reason, we have located the onset point for J/huj = 9 by 
a straight-line extrapolation of the trends in the vicin- 
ity of the half-maximum of the edge curvature feature, 
ignoring the distorted data near the curvature reversal 
itself. 

On the plot of the effective mass as shown in Fig- 
ure ||, the collection of onset points falls into a remarkably 
straight line begging for extrapolation to small values of 
J/Ktv. The fit line describing these points is within error 
of tangency with the J/hui — 1/4 effective mass curve at 
g 0.43, which presents a challenge for interpretation. 
The line of onset points cannot cross the J/Tiuo = 1/4 
effective mass curve since there is no onset phenomenon 
for J/fouj < 1/4. Yet, a true tangency of the onset line 
with the J/hu — 1/4 effective mass curve would imply 
the existence of a curvature reversal for J/hui — 1/4 near 
g w 0.43 that is not evident in our data (see Figure ||). 
The absence of such evidence suggests that within a very 
small interval of J/hw around 1/4, g^ vanishes with a 
very strong, if not singular, dependence on | J /Tito — 1/4|. 

Indeed, the band edge structure is quite singular at 
weak coupling. It can be shown exactly that when 



0+ 



= -2 J J/hu<l/4, (17) 

9=0+ 

= J/hw > 1/4 , (18) 



revealing a true jump discontinuity in the edge curva- 
ture at the crossover from the non-adiabatic to the adi- 
abatic regimes. This discontinuity can be approached 
from J/hui < 1/4 using WCPT, which shows that the 
leading correction to the band edge curvature diverges 
as J/hw — > 1/4 (see Figure |J). This arbitrarily strong 
flattening of the band edge is consistent with the effects 
of quantum mechanical level repulsion as the energy gap 
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between the upper edge of the free exciton band and the 
lower edge of the one phonon continuum vanishes. 




FIG. 6. The breakdown of perturbation theory near 
J/hu = 1/4. Solid lines show edge curvatures according 
to WCPT as in (Q, dashed lines according to SCPT as 
in (^). The solid arrow indicates the limiting dependence 
—2J/huj + ag 2 with a — > oo. SCPT is quite accurate for all 
g at J/hui = 0.1; however, both WCPT and SCPT deviate 
significantly from correct behavior at J/hui = 0.25. Unit of 
energy is Hoj. 

We thus infer that 

9N\j/n u=1 /t = 0; (19) 

however, as a value inferred and not directly computed, 
we distinguish this point on plots throughout this paper 
by a solid square symbol. 



B. Bandwidth 



The polaron effective mass is a pure zone-center quan- 
tity, always characteristic of inner-zone (\k\ < k c ) struc- 
ture; the band edge curvature is a pure zone-edge quan- 
tity, always characteristic of outer-zone (|«| > K c ) struc- 
ture. Although these quantities certainly provide crucial 
information on polaron band shape, and although that 
information is rendered in accurate, quantitative terms, 
it is information that is of a qualitative nature when it 
comes to characterizing the whole band. The polaron 
bandwidth [E(w) — E(0)] on the other hand, relates in- 
ner and outer zone structures and provides an important 
quantitative linkage between the self-trapping transition 
as here quantified in the polaron effective mass, and the 
onset of band narrowing as here quantified in the band 
edge curvature. (See Figure Q.) 
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FIG. 7. Polaron bandwidths {E(ir) - E(Q)]/hu> (32 sites) 
as computed by the Global-Local method vs. exciton-phonon 
coupling strength g; curves left to right: J/fuj = 0.1, 0.25, 
0.5, 1.0, 2.0, 3.0, 4.0, 6.0 7.0, 9.0. As elsewhere in this pa- 
per, square symbols indicate yjv values and diamond symbols 
indicate gsr values. The dashed line at upper left is neither 
derived nor computed, but represents an extrapolation of our 
J/Tioj = 1/4 data to the weak coupling limit. 

The trend evident in the bandwidth sequence illus- 
trated in Figure is typical of the general case; energy 
bands of the adiabatic regime are "clipped" over a non- 
trivial interval of coupling strength, after which a smooth 
but fairly rapid narrowing process commences, culminat- 
ing in the self-trapping transition. From this perspective, 
the self-trapping transition appears to frame one side of a 
distinctly two-sided picture. It is evident that the "knee" 
coinciding with the self-trapping transition is paired with 
another knee at finite coupling coinciding with the onset 
of band narrowing. 

It is clear from the trend in the onset points included 
in the upper portion of Figure [7] that the value of the po- 
laron band width at the onset of band narrowing tends to 
unity as J /Tito — ► 00. Similarly, it is clear from the trend 
in self-trapping points in the lower portion of Figure |?] 
that the value of the polaron band width at self-trapping 
tends to zero as J/htu — > 00. It is evident, therefore, that 
the change in the polaron band width from the onset of 
band narrowing to the completion of this process at the 
point of self-trapping is of order unity in the adiabatic 
regime, and equal to unity in the adiabatic limit. This 
fact will be of particular significance in the following sec- 
tion. 



VI. PHASE DIAGRAM 

The self trapping line and the onset of narrowing line 
divide the polaron parameter space into distinct regions, 
effectively constituting phase boundaries on a polaron 
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phase diagram. For the purposes of this section, we are 
interested in boundaries that transect the whole of the 
polaron parameter space in a manner that a finite data 
sample cannot; thus, we turn our attention to empirical 
curves that accurately reflect our numerical data, but 
permit the consideration of trends beyond the reach of 
direct numerical study. 

Our effective mass analysis of Section IV yields self- 
trapping points that conform well to the empirical self- 
trapping line set forth in previous work [^,^) , namely 

9 st = 1 + y/J/hw (20) 

(see Figure ||). 
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FIG. 8. Phase diagram. Left: Polaron phase diagram according to the Global-Local method, presented directly in terms of 
the "non-adiabatic" parameters J/Tkv and g. Right: Polaron phase diagram according to the Global-Local method, presented 
in terms of the "adiabatic" parameters 7 = Twj/2J and A = g 2 hu/2J. Upper curves (dashed, left (^o|), right (p^)) result from 
analyzing the polaron effective mass. Lower curves (chain-dotted, left (|2l[), right (pq)) result from analyzing polaron energy 
band at the Brillouin zone edge. 



The onset of narrowing line, on the other hand, has 
never before been considered as a construct, so it is nec- 
essary to develop an empirical curve consistent with our 
data. This can be done by considering several character- 
istics of the onset data: i) there is a high degree of par- 
allellism between the onset data and the empirical gsr, 
suggesting that the leading behavior of the onset curve 
may be given by gsr itself, ii) the onset data appears 
to trend weakly toward gsT with increasing adiabaticity. 
iii) the onset curve rises sharply from zero at J/Ulo = 1/4. 
These characteristics are quite constraining, such that by 
considering a construction in simple differences and pow- 
ers we have been led to the functional form 

1 

8 

(21) 

for J/fiuj > 1/4 (see Figure ||), in which there is little 
flexability in the parameters. This function was not ob- 
tained from any formal analysis, and no fit was performed 



On 



= 1 + yJTJfwj- 



8 ( J/Twj - - 
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to determine the parameters in it beyond observing that 
simple whole numbers appeared adequate to provide a 
very satisfactory description of our onset data. 

In terms of the adiabatic scaling scheme, the self- 
trapping curve (p0|) and onset curve (|2l|) take the form 



the latter for 7 < 2, and are shown in Figure Trends 
in the two forms can be compared by noting that in both 
diagrams, fixed- J /Two transits are vertical lines, and cou- 
pling increases from bottom to top. 

The single greatest message of this phase diagram is 
that there are not merely two clearly distinguishable po- 
laron regimes, i.e., large and small, but three. The strong 
coupling region (g > gsr) coincides unambiguously with 
the usual notion of small polarons as highly localized 
quasiparticles with very narrow, nearly cosine-shaped en- 
ergy bands. In conventional usage, the region below the 
self-trapping line {g < gsr) would be termed alternately 
either the "large polaron" regime, "free" regime, "quasi- 
free" regime, or something similar embracing the entire 
interval < g < gsT- Our division of this regime into 
two physically distinct regions, namely, an intermediate- 
coupling regime (<?n < g < gsr) and a weak-coupling 
regime (g < Qn), involves drawing distinctions that run 
contrary to some widely-held views, particularly those of 
adiabatic theory. 

It is perhaps easily enough accepted that the weak- 
coupling regime (g < gjy) should be characterizable as 
the "large polaron" or "quasi- free" regime. What is less 
clear (indicated preliminarily in Ref . |57[ , and expanded 
upon in Refs. Q| and |6|) is that the polaron structures 
found in this regime do not correspond to the large po- 
laron structures familiar from adiabatic theory in 1-D. 
Succinctly put, when g < g^, the characteristic length 
scale for polaron structure is not a ratio such as the 
2J/g 2 h~Lu common in adiabatic theory, but a coupling- 
independent scale such as \J J/Hlu. Thus, the "large po- 
laron" regime we delineate is not occupied by the large 
polarons in the sense of adiabatic theory in 1-D, but by 
states more rigid against the effects of electron-phonon 
interaction. 

The intermediate regime is also not occupied by large 
polarons in the adiabatic sense. Polaron states in the in- 
termediate regime are transitional structures, reflecting 
the realignment of internal correlation structure from the 
stiffcr, weak-coupling structure more characterized by a 
coupling-independent length scale to a softer, strong- 
coupling structure more consistent with the expecta- 
tions of adiabatic theory. This transitional nature of the 
intermediate-coupling regime is clarified by considering 
the adiabatic limit. 




(22) 



(23) 
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A. Adiabatic limit 



While we have stressed throughout that the physical 
phenomenon we call the self-trapping transition should 
be smooth at finite J/Huj and g, it is clear in essentially 
all our data that polaron properties undergo increasingly 
rapid change in the vicinity of this transition as J/hu) and 
g grow large. Just how this steepening trend manifests 
itself depends on the property considered and the manner 
in which the adiabatic limit is approached. 

In the non-adiabatic scaling scheme, the phonon fre- 
quency Hoj is conveniently regarded as the fixed unit of 
energy and the adiabatic limit approached by allowing 
J/hu and g to diverge in the asymptotically-fixed rela- 
tionship g 2 ~ J/huj. In the adiabatic scaling scheme, 
the transfer integral J is regarded as the fixed unit of 
energy and the adiabatic limit approached by allowing 
the phonon frequency to vanish (7 — > 0) as the effective 
coupling parameter A = g 2 huj/2J is held fixed. In the 
latter terms, the self-trapping transition in the adiabatic 
limit is associated with a critical point at 

Ac = \ , (24) 

which can be understood in the non-adiabatic scaling 
scheme as identifying the constant of proportionality be- 
tween diverging g 2 and J/hu> to be unity, consistent with 

A convenient property with which to illustrate behav- 
ior around the adiabatic critical point is the polaron band 
width. As discussed in Section V, with increasing adia- 
baticity, the narrowing of the polaron band occurs essen- 
tially completely within a narrowing interval of coupling 
strength Ag between the onset of band narrowing and 
the completion of this process at the self-trapping tran- 
sition. The steepness of this change can be viewed as the 
ratio of the change in band width AB to the change in 
coupling strength through this process; i.e., 

ABIhuj 1 , s i 

— -f (8J)» -► 00 . (25) 

Ag gsr - gN 

Thus, in terms of non-adiabatic scaling, the steepness of 
the narrowing transition diverges in the adiabatic limit. 
The adiabatic scaling scheme uses a different unit of en- 
ergy that obscures this effect, revealing no divergence; 



AB/J 2 7 3 3 

— t^-~t L_-~ 23 7 s->0. (26) 

This shows that the notion of self-trapping as a criti- 
cal event is not unfounded, but is a limiting aspect of a 
behavior that away from the critical point (finite J/Tiuj, 
nonzero 7) is not associated with a discrete transition, 
nor even with a smooth transition along a discrete self- 
trapping line as described by gsr- Rather, the picture of 
the self-trapping transition that emerges from this work 
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is that of a broad transition involving the whole polaron 
band, having a definable beginning and definable end, 
finitely separated at finite parameter values. Our con- 
structions identify that finite transition interval as the 
"intermediate" regime included between the curves gjy 
(Ajv) and g S r {Xst)- 

We note for emphasis that the entire region included 
between the onset of narrowing (gw or A at) and the self- 
trapping line (gsT or Xst) maps into the adiabatic crit- 
ical point in the adiabatic limit. The polaron structures 
within this region are thus truly "transitional" struc- 
tures, more associated with the adiabatic critical point 
than with any "strong" or "weak" coupling structures 
found in the adiabatic limit. It would appear, therefore, 
that interpretion of these intermediate polaron structures 
in terms of adiabatic theory might fail to faithfully cap- 
ture their transitional character by in effect extrapolating 
from limiting behavior. Similarly, it would appear that 
the nature of the self-trapping transition itself is not fully 
captured in the adiabatic approximation. 

We can make no claim that our curves g$r and <?at 
(Xst an d X N ) are unique in demarking a finite transi- 
tion region that maps into the adiabatic critical point in 
the adiabatic limit; different rapidity criteria will pro- 
duce slightly different traces conveying essentially the 
same physical implications. What is significant about 
the boundaries we have identified is that they are deter- 
mined by physically-meaningful criteria that are model- 
independent, and in the manner of their placement, at 
least, independent of any particular peculiarities of the 
theoretical methods employed in generating the under- 
lying energy band data. That such model-independent 
criteria, applied to high-quality but necessarily inexact 
data should yield a picture of such high internal and ex- 
ternal consistency is compelling. 

B. Adiabaticity 

Throughout this paper, we have used the term non- 
adiabatic to mean J/Tilo < 1 /4, adiabatic to mean J /Tiw > 
1/4, and adiabatic limit to refer to an extreme limit 
J/Tilu — > oo. Our data shows, however, that in many 
aspects there appears to be a crossover from a weakly 
adiabatic regime to a strongly adiabatic regime such that 
it is well to consider the question, "How adiabatic is adi- 
abatic?" . 

We regard the strongly adiabatic regime as that over 
which the trends in all polaron characteristics are quali- 
tatively consistent with the asymptotic behaviors found 
in the immediate neighborhood of the adiabatic limit. 
With decreasing adiabaticity, however, there comes a 
point when these asymptotic behaviors become suffi- 
ciently augmented by contributions ultimately connected 
with non-adiabatic behaviors that the trends in polaron 
characteristics deviate qualitatively from those found 
near the adiabatic limit. This weakly adiabatic regime is 
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not necessarily associated with any limiting property or 
approximation. 

As a broad, smooth crossover, the division between 
the weakly and strongly adiabatic regimes is not sharply 
definable, but we can roughly locate this crossover in 
several ways: 

i) The possible disappearance of an inflection feature 
in the effective mass between J/Tiuo = 1.0 and 1.5 (see 
Figure ||) . 

ii) The extremum in the bandwidth change AB be- 
tween <?jv and gsr between J/Tiuj — 1.0 and 2.0 (see Fig- 
ure @). 

iii) The extremum in the locus of band edge curvatures 
between J /two = 1.0 and 2.0 (see Figure ||). 

iv) The extremum in A at at 7 w 0.38, J/fku « 1.31 (see 
Figure ||b) . 

The aggregate of all these observations would appear 
to imply that the trends consistent with the asymptotic 
behavior associated with the adiabatic limit are not well 
established until J/Tuj > 2 (7 < 1/4), even in those cou- 
pling regimes where " adiabatic" results can be expected 
to hold. Evidently when J/Tilo < 2, or when g < g^, one 
should be wary of applying concepts, approximations, or 
methods tied to the adiabatic approximation. 



VII. DIMENSIONALITY 



The intrinsic complexity of the polaron problem in 1- 
D is not reduced by increasing the real space dimension 
to two or three. Indeed, well-known and widely-invoked 
results tied to the adiabatic approximation suggest that 
polaron structure in 2-D and 3-D should be qualitatively 
distinct from that found in 1-D, and that even the notion 
of self-trapping as discussed in this paper should take on a 



different meaning in higher dimensions 1 34 - 49[| . The root 



of this lies in stability arguments suggesting that in 1-D 
all polaron states should be characterized by finite radii, 
while in 2-D and 3-D polaron states may have either infi- 
nite radii (" free" ) or finite radii (" self-trapped" ) , with the 
self-trapping transition taken to mean the abrupt transi- 
tion from "free" to "self-trapped" states. 

This adiabatic perspective on polaron structure is not 
fully supported by our direct variational calculations even 
in 1-D, strongly signalling a fundamental difficulty with 
adiabatic arguments. Although reportable results from 
our own calculations in higher dimensions are not yet 
available, our 1-D results are quantitatively consistent 
with a variety of independent approaches that also do 
not fully support the adiabatic perspective in higher di- 
mensions. 

Although we must proceed here in a more speculative 
vein than in the balance of this paper, we can reasonably 
infer certain characteristics of polaron structure in 2-D 
and 3-D by interpreting our 1-D results together with 
independently ascertainable dimensional relationships. 
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Using weak coupling perturbation theory, one can show 
that the form of the (isotropic) polaron energy band at 
weak coupling in D dimensions is given by 

E(k) = E^ (k) + E^ (k) + 0{.g 4 } , (27) 

where 

D 

£?(°)(£) =^[-2Jcos(A; i )] , (28) 
i=i 

/•oo 

E W(n) = -g 2 h 2 to 2 / dt e [s (0 '(«)-M« [I (2Jt)f . (29) 
Jo 

From this, it is easily shown that the reciprocal effective 
mass in any direction is given by 

/■oo 

^ = l-g^co 2 J dtte- h ^[e- 2Jt I (2Jt)] D , (30) 

where mo is the free electron mass in the same direction. 
(The generalization to the anisotropic case is straightfor- 
ward.) Comparisons of these weak-coupling results with 
our 1-D masses are shown in Figure 0. 
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FIG. 9. Polaron effective mass, inverse ratio mo /m* . Solid 
curve: Global-Local mass for D — 1. Chain-dashed curves: 
WCPT masses for D = 1, 2, 3, and oo. Scatter-plot: QMC 
data for D — 1, 2, and 3. QMC data kindly provided by P. 
E. Kornilovitch (p|[l| . 

These relations show that at leading order in g, the 
(finite) dimensionality dependence of the polaron bind- 
ing energy (~E<-V(0)) and the effective mass correction 
((mo /m*)— 1) can be quite straightforwardly understood. 
There is no ambiguity, for example, that non-trivially 
correlated polaron states have lower energy than free 
electron states in every dimension, and that these states 
are characterized by effective masses heavier than the 
free electron mass. Elsewhere, we show that these non- 
trivially correlated polaron states of the weak-coupling 
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regime are further characterized by finite radii in any di- 
mension. 

Weak-coupling perturbation theory thus presents a 
picture of polaron structure that is qualitatively similar 
in 1-D, 2-D, and 3-D; in particular, there is no qual- 
itative distinction evident between polaron characteris- 
tics in the 1-D case and those found in higher dimen- 
sions. This qualitative non-dependence on dimensional- 
ity is not limited to the leading orders of weak-coupling 
perturbation theory, but is also evident in the results of 
non-perturbative methods such as quantum Monte Carlo. 
It has been found, for example, that over large ranges 
of electron-phonon coupling strength (including the self- 
trapping transition) quantities such as the electron ki- 
netic energy, the spatial extent of electron-phonon corre- 
lations, and the polaron effective mass behave in similar 
fashion in 1-D, 2-D, and 3-D. In particular, transition be- 
havior that can be associated with self-trapping appears 
to occur in qualitatively similar fashion in each dimen- 
sion, and at comparable values of the electron-phonon 
coupling strength. In detailed comparisons of our results 
with those of other high-quality non-perturbative meth- 
ods, we have found broad and deep agreement; thus, we 
expect that calculations by our methods in higher dimen- 
sions (in progress ||) will produce results qualitatively 
similar to those found in this paper in 1-D. 

There are more intrinsic reasons, however, to expect 
a qualitative dimension-independence of polaron struc- 
ture. We have demonstrated here in 1-D that over the 
entire weak-coupling regime, defined now as the range 
< g < §n, overall polaron band structure is strongly 
characterized by the wave vector k c . This wave vector 
is particularly deeply connected with the value of the 
polaron effective mass since k c persists as a factor in de- 
termining inner-zone band structure to significant cou- 
pling strengths. With k c is associated a length scale kJ 1 
that one might expect to be manifested in polaron struc- 
ture. Moreover, the presence of this length scale would 
likely be more obvious at larger values of J/hu where 
the clipped nature of the weak-coupling energy band is 
more pronounced; in such strongly adiabatic regimes, 
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y/ J/huj. Such a length scale differs markedly 
from that expected under the adiabatic approximation, 
where the polaron radius is generally expected to scale as 
A -1 ~ 2J/g 2 huj. As has been pointed out in a in limited 
fashion in Ref. (57] and more extensively in Refs. [Q and 
H , the length scale characterizing weak-coupling polaron 
structure at large J/Tllo in all our 1-D variational calcula- 
tions is more characteristic of n~ x than of the adiabatic 
scale A -1 . This suggests strongly that the adiabatic ap- 
proximation breaks down at some value of g greater than 
gN, even in 1-D. 

Though we do not yet have comparable variational 
data in 2-D or 3-D, the foregoing arguments general- 
ize. In any dimension, we expect weak-coupling polaron 
structure to be strongly influenced by length scales like 
k" 1 associated with the penetration of the free-electron 
energy band into the phonon continuum. As in 1-D, 
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we expect such weak-coupling structure to persist sig- 
nificantly in any dimension from g — up to a non- 
trivially finite coupling strength g^ associated with the 
the onset of band narrowing and ultimately with the 
self-trapping transition. The polaron structures in this 
regime are not determined by adiabatic stability argu- 
ments, and the length scales characteristic of them differ 
accordingly. Self-trapping is not determined by adiabatic 
stability, but involves the emergence of adiabatic stability 
as a significant factor in determining polaron structure. 

VIII. CONCLUSION 

In this paper we have presented high-precision numeri- 
cal results for the polaron effective mass spanning orders 
of magnitude in both adiabaticity and electron-phonon 
coupling strength, drawn from nearly 1200 complete po- 
laron band structures determined by the Global-Local 
variational method. In addition, selected complete en- 
ergy band structures and other measures of overall po- 
laron band distortion have been examined in detail to re- 
veal significant structure not evident in the effective mass 
alone, nor in the ground state energy or other zone-center 
properties. Resolution of the critical features of polaron 
band structure in their dependence on adiabaticity and 
electron-phonon coupling strength has been sufficient to 
set significant constraints on the interpretation of these 
results. 

Using the effective mass as the primary observable po- 
laron property, we have posited an objective and method- 
and model-independent criterion for locating the self- 
trapping transition based on the rapidity of the rise in 
the effective mass in the intermediate coupling regime. 
The self-trapping line so determined agrees well with that 
obtained by applying similar rapidity criteria to other 
polaron properties, and with a simple empirical curve, 
9st{J /nui), that appears to accurately locate the tran- 
sition regardless of which polaron property is used as a 
favored criterion 

Similar criteria applied at the Brillouin zone edge 
rather than at the zone center provide complementary 
information regarding the onset of polaron band narrow- 
ing. An onset line, g n {J / TUo) , marks the beginning of a 
broad transition phenomenon for which the traditionally- 
understood self-trapping transition at gsr marks the end. 
We are thus led to reinterpret the notion of self-trapping 
as a whole-band phenomenon of finite breadth, spanning 
the coupling interval {9n,9st)- 

The self-trapping line and the onset line divide the 
polaron parameter space into three regimes characteri- 
zable as a small polaron regime {g > gsr), an inter- 
mediate regime (g^ < g < gsr), and a large polaron 
regime (g < g_/v)- Of these three regimes, however, 
only the small polaron regime is unambiguously consis- 
tent with the characteristics traditionally associated with 
these common terms. The states of the intermediate 



25 



regime are transitional structures that map into the adi- 
abatic critical point in the adiabatic limit, and the states 
of the large polaron regime exhibit structural properties 
inconsistent with well-known expectations of adiabatic 
theory, even in the adiabatic limit. 

These characteristics of our findings appear to gener- 
alize to higher dimensions. 

These findings collectively point to a need to creatively 
reconsider the nature of large polarons in 1-D, 2-D, and 3- 
D by means suited to the regime below the onset of band 
narrowing (J/hu large and/or g small, such that g < <7jv)- 
Some well-known methodologies are excluded from this 
regime, including non-adiabatic theories, much of adi- 
abatic theory, most localized-state variational methods, 
and strong-coupling perturbation theory. In work to be 
presented elsewhere, we undertake such reconsideration 
by means of weak-coupling perturbation theory f30|l and 
band-theoretic variational methods such as that used 
in this paper. By directly probing the spatial structure of 
electron-phonon correlations, properties such as the po- 
laron radius can be given precise definition and, in com- 
bination with other polaron properties, used in a con- 
trolled and interpretable way to distinctly characterize 
both large and small polaron structures. 
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